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Abstract 

The small energy anomaly in the single particle density of states of disordered 
interacting systems is studied for the zero dimensional case. This anomaly 
interpolates between the non-perturbative Coulomb blockade and the pertur- 
bative limit, the latter being an extension of the Altshuler-Aronov zero bias 
anomaly at d = 0. Coupling of the zero dimensional system to a dissipative 
environment leads to effective screening of the interaction and a modification 
of the density of states. 
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I. INTRODUCTION 



The density of states (DOS) of disordered interacting electronic systems exhibits singu- 
larity near the Fermi energy. This phenomenon manifests itself in the suppression of the 
tunneling conductance at small external bias and is referred to as zero bias anomaly (ZBA). 
It has been explained and investigated theoretically by Altshuler and Aronov (AA) in 1979 
for d = 3 and latter extended for lower dimensionality by Altshuler Aronov and Lee [Q](for 
a review see Ref. 0). They have predicted that for a d-dimensional system with diffusive 
disorder and weak short range interactions the singular correction to the DOS behaves as 
Sud ~ (max{e, T}) d / 2 ~ 1 with the dimensionality d = 1,2,3. For three dimensional (d = 3) 
samples the ZBA constitutes a small but nonanalytic correction to the DOS (as long as the 
dimensionless conductance, g, is large, g ^> 1). It may be argued subsequently that this 
phenomenon is a precursor of the "Coulomb gap" , which appears in the Anderson insulator 
regime, g < 1. Throughout this paper we shall consider good conductors, hence we shall 
not discuss here the localized regime. But remarkably enough, even for good metals, g ^> 1, 
there is a situation where the ZBA is a strong effect. This is the case with a finite size 
system. The original treatment of AA was given for infinite (quasi) d-dimensional systems. 
For a finite size system the quantization of the spectrum of the diffusive modes becomes of 
crucial importance. For instance, the electron-electron interaction term in the Hamiltonian 
takes the form 

H mt = W V (Q)--PQP-Q--> (!) 

A Q 

where Q = 27rn/L are quantized momenta, pq is the Q component of the electron density 
operator and : . . . : denotes normal ordering. The quantization of the momentum has two 
major consequences for the ZBA. 

First, if one accounts for all finite Q, excluding the Q = mode, in the interaction 
Hamiltonian, the singularity in the DOS is rounded off on the scale max{e, T} E c , where 
E c is the Thouless correlation energy. As long as the ZBA is still a small correction to 
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the total DOS on the scale E c , 6vd(E c ) <C the finite Q modes do not enhance the 
singularity for e, T < E c . It is easy to see || that the requirement 5vd(E c ) <C is satisfied 
for g ^> 1; in this case the finite Q contribution to the ZBA is small for any energy and 
temperature. It means that for such systems the finite Q interaction may indeed be treated 
by perturbation theory, which produces regular expansion in powers of 1/g. Throughout 
this paper we shall assume that the condition g ^> 1 is satisfied, hence the perturbative 
treatment of the non-zero modes is applicable. 

The second consequence of the spectrum quantization is the special role played by the 
Q = term in the interaction Hamiltonian. Its contribution to the energy may be written 
as V(0) : N 2 : /2, where N = p is the total number of electrons in a dot. This term in 
fact corresponds to the classical charging energy of the dot. It leads to a strong singularity 
in the DOS. The first order in interaction perturbative result is 5i>o(e = 0) = — V(0)/(4T) 
(for the Q = contribution T and e are not interchangeable). In the limit T < V(0), the 
perturbative treatment is not sufficient. We present here a relatively simple method which 
allows us to treat Q = contribution to the ZBA exactly. An exact solution is possible 
due to the fact that the zero-mode interaction term commutes with the total Hamiltonian. 
As a result the problem is trivial, although the zero-mode interaction has some interesting 
consequences. The thermodynamical and the response functions are practically unaffected 
by the zero-mode interactions (apart from modifying the statistical ensemble from grand 
canonical to canonical). The single particle DOS is modified in an essential way. We use 
imaginary time functional integral to integrate out all fermionic degrees of freedom. Similar 
methods were used in the context of Josephson junctions with dissipative environment (for 
a review see ||). Non-perturbative treatment reveals the exponential suppression of the 
DOS at the Fermi energy. This result reproduces the classical treatment of the "orthodox" 
theory of the Coulomb blockade 0. In other words, the d — generalization of the AA 
ZBA, and the Coulomb blockade, are two limiting cases of the same theory. 

This statement is in fact not new. In his original paper, Ref . || , Nazarov had established 
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the connection between the two (including a non-perturbative treatment of the finite Q 
modes). A similar, although more transparent, theory has recently been put forward by 
Levitov and Shytov 0. To a large extent, our results can be extracted as the d = limit 
of the expressions found in Refs. The point is that Refs. [|,|| used some uncontrolled 

(although plausible) approximations. We have shown that for Q = all calculations may 
be done exactly, confirming some of the approximations made in Refs. [f)]||. 

One may combine the zero and finite Q mode interactions to the ZBA. In doing so we 
take into account the zero mode in an exact fashion and add on top of it the perturbative 
contribution of the Q / modes. As was mentioned above the latter can be expressed as a 
regular expansion in powers oil/g valid at any temperature energy and dimensionality. To 
address the insulating regime one should consult Refs. 0,|J. In the metallic limit the Q ^ 
contribution modifies the DOS at high energies (e > e 2 /2C), as a precursor to the Coulomb 
blockade which dominates at lower energies (cf. Fig. ||b). 

The new ingredient in our analysis is the inclusion of the zero-point (and thermal) motion 
of the electromagnetic environment in contact with the system. We show that the influence 
of the environment may lead to a Debye screening of the zero-mode interaction, hence to 
a softening of the Coulomb blockade. The physics behind the zero-mode screening is the 
following: the total charge on the dot, eN, interacts with charges in the environment, leading 
to their redistribution. The polarized charge of the environment reduces the energy cost of 
adding (or removing) an electron to (or from) the dot. There is a certain finite time constant 
(the RC - constant of the circuit), characterizing this redistribution of the environmental 
charge, rendering the effective interaction in the dot non-instantaneous (retarded). It might 
appear that the results obtained through this analysis are highly non-universal and depend 
on the particular choice of the model for the environment. We stress, however, that our 
result expressed by Eqs. fl20|) and ([$3|) is quite general, the only model-dependent feature 
is the concrete form of the screened zero-mode interaction, V(u). Nonuniversality in this 
sense is unavoidable, due to the long range nature of the Coulomb interaction: the behavior 
of the dot depends on numerous long distance features. 
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A particularly interesting case is when the zero-mode interaction is fully screened in the 
long time limit. By this we refer to a situation when on a long time scale the addition of 
an electron to the dot does not cost any energy (for instance there is a slow continuous 
leakage of any extra charge from the dot). In this case the usual Coulomb gap in the 
spectrum is absent and one obtains the power-law ZBA, very similar to the one known 
from the physics of Luttinger liquid. The exponent is determined by the time scale of the 
environment polarization. This result for a quantum dot coincide with those obtained in 
Refs. |7]|| for a single tunnel junction connected to a linear RCL circuit ||. For a quantum 
dot one may consider the setup with only partially screened interaction, where the gap in 
the DOS at low energies crosses over to a power-law ZBA at larger energies. 

Let us list several aspects of the problem, which we do not consider here. Although 
we deal with finite size systems, we do not consider effects related to the discreteness of 
the single electron spectrum. It means that the mean level spacing, A, is assumed to be 
the smallest energy scale in our problem. We restrict ourselves to the case of good metals, 
where g = E c /A 3> 1. Thus the energy interval of interest, A < e < E c , is wide. Since 
we are not interested in the single electron spectrum quantization, one may choose any 
boundary conditions for the electron wave functions. We prefer to use periodic boundary 
conditions and employ the momentum representation. We also stress that our analysis 
omits the underlaying periodicity of the problem as function of the charge of the positive 
background (or of a gate voltage) which could be manifested as a sequence of Coulomb 
blockade resonances. Instead we restrict ourselves to the case where the background is such 
that the dot prefers to have an integer charge (half way between resonances). One can easily 
generalize the same treatment for any background charge, noting that in the close vicinity 
of a half-integer (near a resonance) a very different treatment, accounting for multiple 
tunneling events [[H]] is required. In fact even far from resonances multiple tunneling may 
be important (e.g "inelastic co-tunneling" [|11[]). For the sake of clarity we restrict ourselves 
to a "golden rule" scenario, where we consider only lowest order processes in the tunneling 
amplitude, postponing the consideration of multiple tunneling to a further publication. 
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The outline of the article is as follows. In Section |TJ we recall the derivation of the ZBA 
given by AA and extend it to d = 0. The non-perturbative treatment of zero-mode interac- 
tion is discussed in Section |TTTj, where we study in some details the case of an instantaneous 
zero-mode interaction and rederive the results of the "orthodox" Coulomb blockade. Section 
P"V| is devoted to the study of screened retarded zero-mode interaction. We show that in one 
particular case the ZBA reproduces the results derived for a single junction connected to a 
linear circuit. 



II. ZERO BIAS ANOMALY 



The derivation of the ZBA in the diffusive interacting systems proceeds as follows @. 
One calculates the single particle (not to be confused with the thermodynamic, dn/dy) 
density of states as function of energy or temperature. The single particle DOS is defined 
as the imaginary part of the trace of the single particle Green function 



vie) — QTrC/(e 

7T 



(2) 



Following AA |2|], we calculate the first order correction in the screened interaction to the 



single particle DOS in a dirty system. The dominant contribution |L2[ comes from the 
exchange diagram depicted in Fig. |I|. 
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FIG. 1. First order interaction correction to the Green function; wavy line - interaction; 
dashed lines - impurity dressing (diffusons). 
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After performing the fast momentum summation one obtains 

(3) 



in) 



where .D is the diffusion constant; 7i n - the inelastic relaxation rate; V(Q, uj m ) - the screened 
Coulomb interaction, given by 

[ViQ,^)]- 1 = [\/[°](Q)]~ 1 + n(Q,^ m ) (4) 

Here V^(Q) is the bare (instantaneous) Coulomb interaction and II(Q,a; m ) is the polariza- 
tion operator of the system. For an isolated diffusive system the polarization is 



z/'°] is the DOS of non-interacting system. 

Let us consider a three dimensional cube with the linear size L with no current flowing 
through the boundaries. The slow momentum Q assumes quantized values 

Q= T n, 

where n is a vector with integer components, including n = (0, 0, 0). If either the energy, e, or 
the temperature, T, is much larger than the Thouless energy, E c = D/L 2 , one may disregard 
the discreteness of Q and perform slow momentum integration instead of summation. We 
shall refer to such three-dimensional, d = 3. Substituting Eqs. (f|) and in 

Eq. (^|) and performing integrations, one readily obtains [fj] (e.g. for T 3> E c ; e) 

where 03 ~ 3.8 • 10 -2 and g = u^DL d ~ 2 1 is the dimensionless conductance of the sample 
(the full T and e dependence may be found in Ref. 0). At small temperature this correction 
exhibits singular (non-analytic) behavior, which is the ZBA. One cannot employ, however, 
Eq. ([]) for e; T < E c . In this case the discreteness of the momentum spectrum begins to 
play a crucial role. Performing momentum summation in Eq. (|3[), excluding the Q = 
contribution, one obtains for <C e; T <C E c 



where a = £ n _^ (27r|n|)- 6 ~ 1.1 . i(r 3 and 

fo(x) = dyy ^2 - tanh^-y^ - tanh^y^ 

with the asymptotic values fo(x C 1) = ir 2 /3 ; fo(x ^> 1) = x 2 . At T E c Eqs. (§) and (^) 
match parametrically. According to Eq. (0) the small temperature singularity in the DOS 
is rounded off due to finite size effects (the finite value of E c ), (Fig. [^). 

Eq. (0) does not account, however, for the Q = contribution to the momentum sum 
in Eq. @. Below we evaluate this contribution and argue that it reflects the physics of 
the Coulomb blockade. The first problem is that the bare Coulomb interaction, V'^(Q) = 
L~ 3 4iTe 2 /Q 2 , is not well-defined for Q = 0. We argue that for finite size samples this 
expression should be regularized at Q ~ 1/L (cf. analogous regularization of the bare 
interaction in the d — 1 case ||), leading to V^(Q = 0) ~ e 2 /L. More precisely we shall 
use 

^ [0] = £ (8) 

where C ~ L is the self capacitance of the sample, which includes some non-universal geo- 
metrical factors (hereafter the interaction potential, V, without momentum index refers to 
Q = 0). In other words, one may regard the finite value of as the result of fast screening 
processes, occurring on the boundary of the sample, which are not included explicitly in 
Eq. @. According to Eqs. (§) and (|5|), the Q = interaction can not be screened, indeed 
n(Q = 0,uj m ) = 0. This amounts to saying that a finite size isolated system cannot screen 
its total charge. This is not the case for a system which is coupled capacitevely to the 
environment. In the latter case the polarization operator should contain an additional term 
arising from the polarization of the environment. This would lead to an effective screening 
of the Q = interaction. We postpone, however, discussion of the screened zero-mode 
interaction till Section and proceed with the interaction given by Eq. (pi). 
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FIG. 2. Perturbative calculation of the DOS. Full line: Q / modes only; dashed line: Q = 
mode is included. Dotted line: AA result without account for finite size effects (d = 3, cf. Eq. (6)). 
The curves are shifted vertically for clarity. 

Substituting now Eq. ([8]) into Eq. (|3p and performing energy summation and analytical 
continuation, one obtains for the zero-mode contribution to the DOS 



2tt 2 T v2 + 2ttT J 



where ^>(x) is the digamma function. For <C T this reduces to 

Sua e 

— prr = cosh — . 

„[o] 4T 2T 



(9) 



(10) 



Note that the dependencies on temperature and energy are very different, whereas for finite 
d energy and temperature play essentially the same role. The temperature (and diffusive 
constant) dependence is compatible with the AA result, bv& ~ j<(d/2-i) ^ extended to 
d = 0. The zero-mode contribution leads to a dramatic singularity ~ T~ l in the single 
particle DOS, Fig. |[ For sufficiently low temperatures, T <C e 2 /2C, this last expression 
cannot be correct. Indeed in this case Eq. (|T0[) predicts 5v ^> z/'°l, a result which is certainly 
nonperturbative. Moreover, the resulting DOS becomes negative. This clearly indicates that 
the first order perturbation theory is not sufficient to treat the ZBA in finite size systems 
at low temperatures. In the next section we shall show how the zero-mode interaction 
can be treated nonperturbativly. Unlike Eq. (|T0"D, the exact result is well behaved at any 
temperature and, in fact, is well known from the "orthodox" theory of the Coulomb blockade 

I- 



III. ZERO MODE INTERACTION 



A. General formulation 

As we have seen in the previous section, the perturbative treatment of the zero-mode 
interactions in a dot at low temperature meets serious difficulties. From another hand, the 
first order (in the screened interaction) result for Q ^ and a not too dirty system, g 3> 1 
(cf. Eq. (D), (0)), is well behaved. Thus we separate out the dangerous contribution of 
the zero-mode interaction and try to treat it non-perturbativly. To this end consider a 
Hamiltonian describing electrons moving in a disordered potential, which interact through 
the zero mode repulsive interaction only 

(11) 

where a+(a a ) is a creation (annihilation) operator of an electron in an exact single particle 
state (which is defined including disorder potential and spin) with an eigenenergy e a ; is 
the bare zero-mode interaction, Eq. (^); N - the charge of the positive background, which 
we assume to be an integer [|13|| . Bellow we shall comment on the case where Hamiltonian 
includes also finite Q interactions, which can be treated perturbatively. 

Before to proceed we would like to stress the following fact. One may argue that 
the interaction term in the Hamiltonian, Eq. (|IT|), is trivial. Indeed it has the form 
H int = V® : N 2 : /2, where N = £ total number of particles in the sys- 

tem. As N commutes with the Hamiltonian, [H , N] = 0, it does not have any dynamics, 
d t N = 0. Thus the interaction term is just a constant added to the Hamiltonian, and seems 
not to have any nontrivial consequences. This is indeed the case when considering thermo- 
dynamical properties or response functions of the isolated dot. It is, however, not the case 
with the single particle Green function, which correlates the amplitude of the creation of 
one additional electron at an initial time, £j, and its subsequent destruction at time £/. Thus 
any measurement of the single particle Green function assumes implicitly the existence of 
processes which do not conserve the number of particles in a dot (e.g. due to tunneling). 
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53 a+a a - N 



Such processes render the entire Hamiltonian noncommuting with the particle number. As a 
result the TV 2 interaction term has a significant impact on the single particle Green function. 
After completing the calculations we shall comment on the relation of the above arguments 
to gauge invariance. 

The imaginary time single particle Green function may be written as [ 13] 

1 



with the fermionic action given by 



(12) 



SW a M = / dr 
Jo 



Z V£ ( T ) ( 9 r + Ca - V)lpa (t) + 



(13) 



here Z{p) is the partition function and \x is the chemical potential. Splitting the interaction 
term in the action by means of the Hubbard-Stratonovich transformation with the auxiliary 
Bose field, <p(r), one obtains 
1 



Gain-Tf) 



( r )i e -/ ^^MpvPii-^W-tiVb^T)] 



(14) 



with the same transformations in Z(ji). Here and Gg\ T ii T fi /-0 are respectively 

the partition and Green functions of non-interacting electrons in the time dependent (but 
spatially uniform) potential, i(p(r). To calculate these quantities one should resolve the spec- 
tral problem for the first order differential operator [d T — fi + e a + i4>(r)] with antiperiodic 
boundary conditions. This can be easily done (see e.g. chapter 7 in Ref. ||14]| ). For the 
spectral determinant (partition function) one finds 



z^) = z^-uM, 



(15) 



where Z^\ji) = exp{— /3f2'°l(//)} is the partition function of non-interacting elec- 
tron gas. We have introduced Matsubara representation for the boson field, 0(r): 
(f) m = f£ d,T<p(r) exp{ia; m r}, u m = 2irmT. The Green function is given by 
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Gf (r h r f , ai) = 0j? (r< - r,, „ - i^tf (16) 



where C^(e n ,/i) = (ze„ — e Q + /i) 1 is the Green function of non-interacting fermions. It is 

convenient to rewrite the exponent in the last equation in the following form 

T n,T f 
i" — 

where Jm' T/ is the Matsubara transform of the following function 



exp{* P dr[<P{r) - O ]} = exp{/3 £ ^Ln^ }} (l 7 ) 



^.'■/ = < y(r-r < )-<J(r-r / ). (18) 

Transforming next the functional integral over 0(r) to integrals over the Matsubara 
components, m , we obtain 

<?«(r, - r f ) = -±-J d^ e -^^ m ]-^^o+n™&^o)]^M( Tl _ T/(/< _ , 0O ) 



y n ^ ex p \ p Yi 



2Vt°l 



(19) 



m^O I m^O 

Again with the analogous modifications in Z{\x). The integral over the static component, 
00; describes the smooth transition between the grand canonical ensemble with the chemical 
potential /x (at V' ' = 0) to the canonical ensemble with No electrons (at = oo). For 
large enough systems (A -1 = — <9 2 f2^/<9/i 2 3> 0) one can neglect differences between the 
two statistical ensembles. This means that the integral over 0o can be calculated in a saddle 
point approximation, leading to { T i ~ T fiV)i where the stationary point, /Z, is the real 
solution of the equation {p, - ]£) /V® + N + dQ^ (p) /<9/i = 0. The remaining integrals (over 
4> m for are purely Gaussian. As a result one obtains 

g a ( Ti - r /; //) = 0M(r< _ TflJ l)e- s ^f\ (20) 

where 

S(t)=TJ2 — (l-e 1 ^). (21) 

Eqs. (^Up and fl2"T| ) solve the problem of finding the exact single particle Green function in 
the presence of the zero-mode interaction. This result is depicted diagrammatically in Fig. 
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|3|. The dressed Green function of the interacting problem is given by the bare one (with a 
renormalized chemical potential) decorated with the propagator of the auxiliary boson field 
T>(t) = exp{— S(t)}. We refer to this auxiliary field as a - Coulomb boson. As T>(0) = 1, 
the zero mode interaction does not influence equal time Green functions (apart from the 
renormalization of the chemical potential). As a result thermodynamical quantities are not 
affected by the presence of the Coulomb boson. This is not unexpected since the interaction 
term commutes with a total Hamiltonian, hence it does not affect quantities defined for the 
closed system. By contrast, to measure a two point Green function one should perform a 
tunneling experiment, where the system can not be considered as completely isolated. In 
this case the presence of the interaction term in Eq. (|TTD (and hence of the Coulomb boson) 
is of crucial importance. 

3){u) 
/ \ 

T\\ . i T 2 

FIG. 3. Diagrammatic representation of the Coulomb boson. 



This observation can be related to gauge invariance |15|| . As we have seen the problem 



with zero mode interaction is essentially reducible to that of an electrons gas in a spatially 
uniform a.c. potential. Such a potential can be always removed from the problem by a 
time-dependent gauge transformation. Thus gauge invariant physical quantities (e.g. ther- 
modynamical quantities) are not affected by the presence of spatially uniform potential, cf. 
Eq. dl5|) . The single particle Green function, being non gauge invariant object is affected. 
By introducing an electron tunneling from an external source we fix the (time dependent) 
phase of the electron wave function in the system. This point has been discussed by Finkel- 



stein ||15[| , who argued that unlike the conductivity or thermodynamical quantities, the non 
gauge invariant single particle DOS may be affected by very small Q terms, leading in d < 2 
to more pronounced singularities. 
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B. The Coulomb boson and DOS 



For the instantaneous (frequency-independent) interaction, Eq. (|8|), the sum in Eq. ( |2T|) 
may be easily performed, yielding 

Sir) = - (|r| - -) . (22, 

Using the Lehmann representation for temperature Green functions [ll]] one may write for 
the propagator of the Coulomb boson 

, , f°° duo' B(uj') . . 

J-oo 2ir iuj m — 00' 

where the spectral function B(u) is defined as B(uj) = —2QV r (uj). Fourier transforming of 
T>(t) = exp{— S(t)} and analytically continuing we obtain (see Appendix |A|) 



B (") = )lyWr(e 2vl0]T -e avPiT j. (24) 

According to Eq. (|20|) (cf. also Fig. ||]) the electron Green function is given by 
Ga(tn) = TY,uj m <3aK e n — w m )D(u m ). Performing the summation by a standard contour in- 
tegration and then analytical continuation (ie n — > e + i8) one obtains for the one particle 
(tunneling) density of states, u(e) = —n^ 1 J2 a ^$&a ( e )' 



vie) = - 



- r — (tanh^ + coth^) B(u)^(e - W ), (25) 
2 J-oo 2tt V 2T 227 v ; v ; v ; 



where u^(e) is the density of states in the absence of the zero-mode interaction and B{u) 
is given by Eq. ([24]). For non-interacting electrons (with or without disorder) i/'°l(e) may be 
well-approximated by a constant value, z/'°l, in this case the integral in the last expression 
may be readily evaluated. In the limit of weak interaction (V^ <C T) one obtains 

HZ = 1 " ^osh"^. (26) 
„[o] AT 2T y ' 

This is the zero bias anomaly of AA extrapolated to d = system, which has been obtained 
from a diagrammatic expansion, Eq. QTUp. Notice that we did not actually use the fact that 
the system is disordered. Note also, that the exact lowest order result, Eq. (|2~6"1) , coincides 
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with the one obtained from the exchange diagram only, Eq. ([101) . The first order Hartree 
term with the Q = interaction leads to redefinition of the chemical potential which has 
been absorbed in the factor, fx. For strong interaction (V^ ^> T) one has (Fig. |4]a) 



! 2ttT 



r[0] 



r[0] e~ — cosh— ; e < W®T < V^, 



1-e 



-(e-y[°])/T. 



(27) 



e > V@ > T. 



The exponential suppression of the tunneling density of states near the Fermi energy is a 
direct manifestation of the Coulomb blockade. The ZBA in d — systems, being treated 
non-perturbativly, leads to results which are well-known from the "orthodox" theory of the 
Coulomb blockade, cf. Ref. H. 




(b) 



FIG. 4. Non perturbative DOS. (a) as function of temperature at the Fermi energy (e = 0); 
only the Q = mode is included, (b) DOS as function of energy at T = 0; finite Q contributions 
are included. 



One may ask how the finite Q interactions modify the standard Coulomb blockade pre- 
dictions. In the case where the finite Q interactions can be treated perturbativly the answer 
can be read off Eq. (^): in Eq. (|25|) one should substitute the perturbative AA expression 
for the DOS [y AA ) for z/M(e). Indeed, one may repeat the derivation given above in the 



15 



presence of the finite Q interactions; the result coincides with Eq. (P0|), where Qi°> includes 
the effects of all but zero-mode interactions. In the limit of zero temperature one obtains 

f 0; lei < VW/2, 

i/(c) = (28) 
I v AA (\e\ - V®/2); \e\ > V^/2. 

This result (Fig. |]b) implies that for voltages larger than the Coulomb blockade gap the 
tunneling current is still somewhat suppressed due to the presence of the finite dimensional 
ZBA. Although the finite Q interactions reduce the jump at |e| = V [0] /2 (T = 0), it do not 
remove the discontinuity. There is, however, important physics which leads to the rounding 
off of the threshold at |e| = e 2 /2C even at zero temperature. This is the screening of the 
zero-mode interaction due to fluctuations in the electromagnetic environment. This issue is 
to be discussed next. 



IV. SCREENED ZERO MODE INTERACTION 

In the preceding section we studied the tunneling DOS of a dot which is perfectly isolated 
(both electrically and electromagnetically) from the outside world. In this case the total 
charge of the dot can not be screened. This is why we have employed a bare (instantaneous) 
zero-mode interaction, Eq. @. This scheme, however, is not very satisfactory. First, 
measurements of some quantities of major interest (such as the single particle DOS studied 
here, or the inelastic broadening to be discussed elsewhere), require by definition direct 
coupling to an external medium. But more importantly, even in the absence of external 
leads the dot interacts (through capacitive coupling) with external gates, conducting layers, 
external charges, all of which will be referred to as the environment. The creation of an 
additional charge in the dot leads to redistribution of charges in the environment, which 
takes finite time. The redistributed environmental charge partially screens the initially 
created charge on the dot, reducing the zero-mode interaction energy. The fact that this 
redistribution is not instantaneous implies that the effective screened interaction is retarded. 
Equivalently one may notice that the dot's capacitance is proportional to the dielectric 
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constant of the surrounding medium. Unless this medium is a perfect vacuum, the dielectric 
constant is a function of frequency. As a result the zero-mode interaction, e 2 /C(u>), is not 
instantaneous, but rather characterized by a finite retardation time scale. 



H2> 




U Z(w) C G C 

FIG. 5. Equivalent circuit of a dot coupled to the enviroment. The grey square represents 
a dot, C being its self capacitance; Cq~ the mutual capacitnce to the environment; Z{u)- the 
linear impedance of the environment. Here U represents the noise voltage. The black arrow is a 
tunneling tip through which particles may be injected to the dot, measuring its DOS. This tip may 
be coupled to the dot through a large impedance barrier. 

Phenomenologically it is convenient to consider an equivalent circuit (Fig. ||), assigning 
the self capacitance, C, (bare interaction) to the dot, which in turn is capacitively coupled 
(through Cq) to the gate electrode. The gate is electrically coupled to the ground through 
the linear impedance Z(uj). The voltage source, U, represents the equilibrium noise voltage 
of the entire circuit. We shall model the dot-environment interaction by assuming that 
the dot is subjected to a time dependent (but spatially uniform) noise potential due to the 
environment 

Hnoise = V(t) a t a a- (29) 
a 

At the end of these calculations all physical quantities should be averaged over realizations 
of the noise. We shall assume further that the noise is Gaussian with zero mean value. The 
averaging procedure thus may be written as 

(. • .U = jfj V[ V (r)]e-U f>^v(r)K-HT-r ')„(.') _ _ _ ^ (30) 
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where Af is a normalization factor and K(r — r') = (v( T )v( T> )) * s the noise correlator, which 
is to be determined employing the fluctuation-dissipation theorem (FDT). An important 
observation is that the partition function of the dot, Z((a), is not affected by the noise term. 
Indeed, we have already seen, Eq. flTSP , that the partition function of the electron gas in 
a spatially uniform a.c. field depends on its mean value only (which is zero in our case). 
This can be understood as a consequence of the gauge invariance of the partition function - 
spatially uniform field may be always removed by a gauge transformation. As a result, the 
noise term in Eq. fll2[) enters in the numerator only. Averaging Eq. flTJD over the Gaussian 
noise, Eq. (|30|) , leads to an effective fermionic action with interaction which is non-local in 
time, 

SintiKM = \ If drdr'Y,r a {r)Mr) (V [0] 5(r - t>) - K{r - 0) (31) 

As a result, one obtains an effective renormalization (screening) of the zero-mode interaction 
potential 

V^^V(u m ) = V^-K(u m ). (32) 

Further calculations follow the same steps outlined in Sec. |TJ. The final result is given again 
by Eq. © where now (cf. Eq. @) 



S(t) = T ]T ^^(1 - e^ r ). (33) 

We shall next employ the FDT to calculate the noise correlator, K(iu). According 
to the FDT the equilibrium noise spectrum (T = 0) of the total noise voltage generated 
by the circuit is (UU) = e 2 iuZ tot , where the total impedance of the equivalent circuit is 
Z to t = {iujC)~ l + {iujCg)^ 1 + Z(u). The corresponding voltage drop on the dot is r\ = 
UiiuiC) -1 1 'Ztoti thus the noise correlator, K = (rjr)), is given by 

r _ , e 2 l/(iujC) 

KM = 77 '\ ■ (34) 

Substituting this expression in Eq. ( |3"2] ) and rewriting it in a finite temperature form one 
obtains 



1 m ' C l + C/C a + \^ m \ZC y ' 

The high frequency (unscreened) value coincides with that of Eq. (§), whereas at low 
frequency the interaction is partially screened and is given by the total capacitance, 
V(0) = e 2 /(C + Cg)- The characteristic crossover frequency is given by the U RC" time 
of the circuit, uo k, (ZC)~ l . As a result the long time behavior of S(t), determined by the 
small frequency asymptotic of the screened interaction, is given by S(t) ~ \r\e 2 /2(C + Cg) 
for ZC < r < (3. 

A case of particular interest is that of a "maximally" screened interaction, when 
V(iu — > 0) = 0. In this case the linear term in S(t) is absent and the action grows at 
most logarithmically at large r. This would immediately imply that instead of exponential 
suppression of the DOS at the Fermi energy, one obtains only a power-law ZBA. According 
to Eq. fl35|) such a case may be realized when Cg — > oo (more precisely ZCg 3> (3). This is 
the case of a dot strongly connected to the environment. For simplicity we restrict ourselves 
to the scenario of a pure ohmic environment, Z = R, in which case the interaction potential 
is given by JT7 



VM = V® 7^~~T^ r> (36) 



where Q = (RC)- 1 and V® is given by Eq. (|8|). As V(u = 0) = the addition or subtraction 
of an electron from the dot costs no Coulomb energy over long time scales. On short time 
scales, however, before the environment adjusts to screen out the added electron, one has to 
pay some energy, Eq. (^). Thus the addition of an electron can be considered as "tunneling" 
under an energy barrier in the time direction. The very same physics has been discussed 
in the context of d — 2 systems in Refs. [18||. The related energy cost on short time 



scales suppresses free particle exchange between the dot and the particle reservoir (although 
Coulomb blockade in its strict sense is absent). This leads to the suppression of the tunneling 
DOS hence to ZBA. Substituting Eq. @ into Eq. (0), one obtains for ft" 1 < r < p p 



S(t) w 2r log(l + fi|r|) + O ((fir)- 1 ) , (37) 
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where f2 = e 7 fi with 7 = 0.577 . . . being the Euler constant and r = / (27rfi) = Re 2 / (2nTi) 
is dimensionless resistance of the environment. Performing the Fourier transform and ana- 
lytical continuation (see Appendix [A]), one obtains for the zero-temperature spectral density 
of the Coulomb boson 



B(u) 



sign uj 2tt 



2r-l 



-\uj/q\ 



(38) 



n r(2r) 

Finally the zero temperature DOS may be found from Eq. (|2~5"D with B(u) given by Eq. (|38"1) . 
Assuming a constant bare DOS (i.e. z/M = const) one obtains 

1 

1 - r ( 2r. 




It 



lei < fi, 



|e| > max{Q,y[°l}, 



(39) 



where T(2r, x) is the incomplete gamma function. For the maximally screened zero-mode 
interaction scenario the ZBA in d = has a power law behavior, rather than the gap 



given by Eq. (|28|) . For high-impedance (slow) environment, r ^> 1, there is a crossover 
to the "orthodox" Coulomb blockade in the interval £1 < e < V^. However for the low- 
impedance, r < 1, environment the power-law ZBA, Eq. (|39|) , directly crosses over to the 
finite dimensional AA result at e ~ Q. 

The straightforward finite temperature generalization of Eq. ([38|) is (see Appendix |A|) 



B(u) 



sinhw/2T 2 f2nT\ 2r - 1 ^ ^ + 



T(2r) V Q 



2nT , 



(40) 



Substituting this result into Eq. 
e = 0; T < ft 



one obtains e.g. for the DOS at the Fermi energy, 



^0 



2r 



a(r) hn 



(41) 



where 



a(r) 



7rr(2r) Jo 



00 , |T(r + zz) 
ax- 



il' (±)2r; r«l, 



cosh 7ra 



2- 2 v^F 



r > 1; 



(42) 
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a(l/2) = 1/tt; ct(1) = 1/8. At T > max{fi, V"' '} the DOS obeys Eq. ©. Screening of the 
zero mode interaction converts the exponential suppression of the DOS (cf. Eq. fl27|) ) into 
the power-law, Eq. (|4l|) . 

Our results in the limit of maximally screened interaction Cq —* oo are compatible with 
the results of Refs. where tunneling through a single tunnel junction coupled to a linear 
impedance was considered. This coincidence is not surprising at all; indeed a dot coupled 
to a circuit through large capacitance, Cq, can be viewed as an island from which charge 
can continuously leak out; this is practically the setup of Refs. 0,f|. Our present analysis 
stresses the physics of a weakly coupled dot (with a finite Cq, Eq. fl35|)) and its relation to 
the ZBA. At T = 0, employing the above given expressions one obtains 



!/„(£) | 0; |6|<e72(C + C G ) 



7/ 



(43) 



where r = t\Cq/{C + C^)] 2 - At finite temperature and e = one has 



exp{-e 2 /8T(C + C G )}; T < e 2 /2(C + C, 



g 



^ « <! (T/f2) 2 ^; e 2 /2(C + C G ) < T < fl, (44) 



1 - e 2 /4TC; max{0, e 2 /2C} < T 

(the above expressions account for the Q = contribution only). We stress that, our results 
are a particular, d=0, case of a general non-perturbative expression for the ZBA obtained by 
Nazarov |J and Levitov and Shytov 0. Our point here is that for d = case all calculations 
can be carried out exactly, avoiding some of the uncontrolled albeit plausible assumptions 
employed in Refs. |5]||. 
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APPENDIX A: 

This Appendix is devoted to the calculation of the spectral density of the Coulomb 
boson, V{t) = exp{— 5(r)}. In the case of unscreened zero mode interaction S(r) is given 
by Eq. fl22|). The Fourier transform has a form 



Jo 



(Al) 



Deforming the contour of integration as it shown on Fig. || and making the obvious redefi- 
nition of variables one obtains 



V(Ur. 



i I dte 
h 



2/3 



-iVl°h/2 _ iV^t/2 



(A2) 



In this expression one can perform analytical continuation, resulting in T> r (oj) = T>(u m — > 
—iuj + 5). The imaginary part of T> r (oj) may be easily evaluated extending the region of 
integration up to — oo. As a result, B{uj) = —2QV R (u) is given by Eq. (p4|) . 




Rer 



FIG. 6. Countour of integration in Eq. (Al). 



Next we discuss the maximally screened zero mode interaction. For < r /3/2 the 
action is given by Eq. (|37|). Periodicity and symmetry requires that for (3/2 t < (3 
action has the form S(r) w 2rlog(l + Cl(/3 — r)). Performing the Fourier transform by 
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deforming the contour of integration and analytical continuation exactly as for the non 
screened interaction, one obtains B(uS) = —F{uj) + F(—u), where 

0(w) 2tt fuj^ 2 ^ 1 



F(lu) = 23? 



Aujt 



dt- 



(A3) 



/o "(i + int) 2r n r(2r)\n, 

This leads to zero temperature expression, Eq. (^). To generalize the above expressions 
to finite temperatures we use the conformal transformation r — > smirTr / (nT) (cf.Ref. p0| ) 
yielding 



F(w)^S 23? 



-nrr „oo 



dte 



iwt ( \ 



2r 



tt 2r Jo VsinhTrTt 

n r(2r) KIT) 



r r + 



(A4) 
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